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Abstract 

Matrix quasi exactly solvable operators are considered and new conditions 
are determined to test whether a matrix differential operator possesses a or 
several finite dimensional invariant vector spaces. New examples of 2 x 2- 
matrix quasi exactly solvable operators are constructed with the emphasis set 
on PT-symmetric Hamiltonians. 



1 Introduction 



Quasi exactly solvable (QES) operators refer to a class of linear operators (typically 
of Schrodinger type) which preserve a finite dimensional subspace of the Hilbert space 
on which they act [1, 2, 3]. Is most of the examples known, QES operators can be 
transformed into operators preserving a space of polynomials of given degree after a 
suitable change of variable and change of functions (also called "gauge transforma- 
tions"). The operators preserving a space of polynomials therefore play an important 
role in the study of QES operators. In the case of scalar equations, the gauge trans- 
formation consists in factorizing the ground state out of the wave function, the change 
of variable can be performed in a straighforward way (at least formally because it 
leads usually to elliptic integrals). 

In the case of coupled equations, where the operators appear in the form of a 
matrix whose components are differential operators, the construction of the gauge 
transformation setting the operator in a form which manifestly preserves a vector 
space where component are polynomials, appear more tricky, see e.g. [4, 5]. 

In the second section of this paper, we establish a set of algebraic conditions 
to test wether n x n matrix-valued operators of a certain type preserve a vector 
space of n-uple of polynomials with component of definite degrees. We work with 
2x2 matrices but the method can be extended to higher dimensions. This new 
method was tested on all QES known matrix equations. In the other sections, we 
take advantage of these conditions and construct several new families of QES systems 
where the emphasis is set on PT-symmetric invariance. This original issue for the 
mathematical framework of quantum mechanics was proposed in [6] and developped 
is several subsequent papers but, to our knowledge, it has not been studied in the 
context of coupled systems of Schrodinger equations. 

In Sect. 3 we propose several matrix extentions of the Razhavi operator. Scalar 
Razhavi-types of potentials were considered recently to produce examples of PT- 
invariant, non-hermitian potentials with real eigenvalues [7, 8]. Here, we develop 
matrix extensions of them both with trigonometric and hyperbolic potentials. 

In Sect. 4 we obtain a matrix generalisation of the QES example of PT-symmetric 
hamiltonian with an anharmonic potential of degree four [6] and reconsidered recently 
[9]. Finally, in Sect. 5, we show how the problem of section 4 can be transformed 
into a system of recurence equations in the spirit of [10]. 

2 Matrix QES Operator 

In this section, we propose a general test to check whether a 2 x 2 matrix differential 
operator H (depending of the variable x and of the derivatives d n /dx n ), preserves 
a vector space whose components are polynomials of suitable degrees in x. We will 
consider a family of operators H which can be decomposed according to 

H = H 1 + H + ..., (1) 
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The diagonal components of Hi are differential operators of degree 1 (i.e acting on 
a generic polynomial of degree n in x they increase the degree by one unit) and the 
off-diagonal elements (#1)12 and (#1)21 are respectively proportional to x s and x 5 , 
with 5 = 0, 1, 2 and 5' = 2 — 5. The operators if , • • • have lower degrees in all their 
components than the corresponding components in H\. The dots in (1) represent 
operators with lower degree ordered according to the same rule. 

Most of the QES matrix operators known [5, 13, 12] can be reduced to the form 
above after a suitable gauge transformation and change of variable. The different 
components of Hi, H contain several constant parameters which label the physical 
coupling constants in case where H is an Hamiltonian operator. 

Now, we try to obtain the conditions of Hi and H (i.e. conditions on the coupling 
constants) such that the operator H possesses a finite dimensional invariant subspace 
of polynomials of the form 

V = span J ( Pn )\ , neU , m = n-5 + l (2) 



q. 



m . 



where p n , q m denote polynomials of degree n, m in the variable x. For a generic vector 
in V of the form 

_ / a x n \ ( aix n - 1 \ , . 

%P ~ {p x n - s+1 ) + \Pix n - s ) + ---> [6) 

where a iy Pi are complex parameters, the vector Hip can be decomposed according to 



H^ = diag(x n+ \x n - s+2 )Mi 



( "o 




+ (diag(x n ,x n - 5+1 )M 1 +diag(a; n ,x n - <5+1 )Mo 
+terms of lower degrees , (4) 

where the constant 2x2 matrices M 1 ,Mi and M can be obtained after a simple 
algebra. 

The necessary conditions for V to contain an invariant vector space of the operator 
H read 

where the second condition has to be fullfilled irrespectively of the values oii,fli. 

The condition {%) implies detMi = and the vector (cto, /3qY to be a zero- 
eigenvalue eigenvector of Mi. This fixes the relative coefficient of the terms of highest 
degree in V (see Eq.(2)). The condition (ii) can be fullfilled only if the following 
conditions hold 



PoJ KPoJ ' 1 V «o J V0 
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where M l means the transpose matrix of M. 

The conditions (i),(ii') allow to reconstruct in a systematic way the invariant vector 
spaces of all QES operators presented e.g. in [12, 5, 13]; in particular the conditions 
on the different parameters and the relevant changes of variable now emerge in terms 
of elementary algebra on matrices. 

In order to illustrate this method, we reconstruct the invariant vector space of the 
QES Hamiltonian [12, 13] 

H(y) = -^2 + M 6 (y), (8) 
where M 6 (y) is a 2 x 2 hermitian matrix of the form 

M 6 (y) = {Aply 6 +8 PlP2 y 4 +(Apl-8m P2 +2(l-2e) P2 )y 2 }-5L 2 +(8 P2 y 2 +A Pl )a 3 -8m P2 K a 1 

(9) 

It is known that after the usual "gauge transformation" of H(y) with a factor 

<P(y) = y e ^p-{jy 4 + Piy 2 } , e = o,i (10) 

and the change of variable x = y 2 , the new operator H(x) is obtained 

H{x) = r\y)H{y) ( t>{y)\y=v- x ( n ) 

For simplicity we assume e = 0,pi = in the following. The operator obtained H(x) 
then reads 

^> = <- 4 ^- 2 l^ + 8 ^ J+(m o~ 2) J+ U)- Sm ^ < 12 > 

with J+(m) = x 2 d x — mx. It can be decomposed along the lines of Eq.(l) : 

H(x) = H 1 + H + H_!, (13) 



with 



* =^( J+( r 2) ;;:,) 

Ho =o, 



d 2 d 

E_ x = ("4a;— - 2— )H 2 - 8mp 2 K a_ (14) 

In this case, {H\)i 2 is a constant(i.e. 5 = 0), while (Hi) 2 i = 0, in addition the 
operatot H is zero. The invariant vector space to be looked for is of the form 

_ / aox™- 1 + a x x m - 2 + ...\ , s 

\ (3 x m + (3ix m - 1 + ... ) 
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The determinant of the matrix Mi is trivially zero and the condition (i) implies 
^ = rriKo- The first conditions (W) is trivial since M = (as a consequence of 
H = 0). Finally, the second condition (W) can be easily checked : 

Hi) =-^(T) = -*^(?) 

-8 P ,a/3„(7°) = -a(£) del 

In the following section, we will present several examples of QES matrix operator 
based on extensions of the scalar Razavi potential. 



3 PT invariant non- hermit ian matrix Hamiltonian 

In [7, 8] PT-invariant models based on the scalar Razhavi potential are analyzed with 
the emphasis set on the reality properties of the spectrum. This can be done partly in 
an analytical way because the potentials considered are QES. The authors considerd 
both, hyperbolic and trigonometric cases, invoking an anti-isospectral transformation 
[11] to relate the spectra of both types. Here we will consider matrix extensions of 
these equations and see that several form of the non diagonal elements Hi 2 and H 2 i 
can lead to QES operators. We will first consider periodic potentials, formulated in 
terms of trigonometric functions. Then an example involving elliptic functions for 
the potentials will be presented. 

3.1 Trigonometric case 

From the unidimensional potential studied in [7], we will build a family PT invariant 
matrix Hamiltonian and use the technique developped the previous section to check 
its quasi exactly solvability. We start from a general Hamiltonian of the form 

H= f-£ + (pcos2x-tM) 2 + A H 12 

V #21 -£z + (pcos2x-tM) 2 + A 

where p is a free real parameter and A, A, M, M are constant to be specified. There 
are several forms of H i2 , H 2 i which lead to QES operators. One can assume M > M 
without loosing generality. The general properties of the diagonal component of H 
and of trigonometric fonctions will reveal that QES operators can be constructed by 
choosing if 12 according to one of the following form 

Hi 2 = C cos2x + D or H 12 = Ccosx or H 12 = C sin x or if 12 = C cos x sin x (18) 

and similar forms respectively for H 2 i with, however, a priori independent coupling 
constants for C and D. 
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In order to reveal the algebraic properties of this family of operators, it is conve- 
nient to perform a first gauge transformation according to 



fr -*coe2x(*- e (l-*)- )rr p 0cos 2 X A 1 - zf 

^ z-\\-z)-<t>) \ z\\-z) 

Hn H 12 
H21 H 22 



(19) 



where z = (cos2x + l)/2. Further choosing the parameter 6, according to — i| 
the components of H are obtained after an algebra: 

d 2 d 
H n = -Az{\ -z)— + 2(2z - 1 - 4(1 - z)e + A<pz)— + p 2 - M 2 + 80e + 2e + 20 + A 

dz 2 dz 

-%ip{z{\ - z)-^j- + e(l -z)-<j>z+ M ~ 1 (2z - 1)) 

H 12 =z i -^l-z)^H 12 , 
H 21 = z^(l-zf^H 21 , 

H 22 = -4s(l + 2(2^ - 1 - 4(1 - z)e + ±fa)^j- + p 2 - M 2 + 80e + 2e + 20 + i 

-8ip(z(l -z)-^ + e(l-z)-<f>z + ^^(2z - 1)) (20) 
dz 4 

and where we have neglected the singular terms of the form 

^—^2e{2e - 1) + -^20(20 - 1) (21) 
z 1 — z 

in Hn (and a similar terms with e — > e , — > in H 22 ) since we assume from now on 
e(2e - 1) = 0(20 - 1) = e(2e - 1) = 0(20 - 1) = (22) 

The different choices for H\ 2 proposed in Eq.(18) now appear to be natural since 
they will automatically lead to a polynomial expressions in z when the choice of the 
parameters e, e, 0, is done according to Eq. (22). In the following, we will analyze 
in details the case H X2 = C sin x cos x. The algebraisation corresponding to the three 
other cases can be done similarly. 

In this case, the possible values for the parameters e, e, 0, allow for four alge- 
braisation, for the wave function ip = (ipi,ip 2 ), namely : 

type i : ip = (p n , sinxcosxg„_i) 
type ii : = {p n -\ sin x cos x, q n ) 

type Hi : ip = {p n sin x, q n cos x) 

type iv : ip = (p n cos x, q n sin x) 

where p n , q n , etc. denote polynomials of degree n in the variable z. 
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Acting on an eigenfunction of type (i), the conditions for algebraic solutions are 
e = = 0, e — (j) — 1/2. The operator H can then be decomposed according to the 
prescription of Sect. 2, leading to : 



H = Hi + H + H_i, (23) 



whith 



d 2 /M „. s d 9X1T (A' Cz 



and 



with 



M - ^ -2 * (26) 



M • .V/ l// . (I - 4n 2 ) + MM = T7r ^ (28) 



A' = A - M 2 + 2ip(M — l) , A' = A + 4 — M 2 + 2ip(M - 3) (27) 

Using the parametrisation corresponding to type (i) for the wave function, we 
can easily obtain the form of the matrices Mi, Mi, M and the conditions on the 
parameters leading to QES operators. In the present case, we got 

CC 

The condition involving M fixes the difference between the constant A, A, namely 

A — A = M 2 — M 2 (29) 

Considering the parametrisation corresponding to type (ii) for the wave function, one 
can easily obtain e = 0= l/2, e = = O and the associated operator H. The 
action of H on an eigenfunction of the type (ii) gives after an algebra the matrices 
Mi, Mi, M which lead to the same QES conditions found in the previous case as 
given by the Eq.(28) and the Eq.(30). 

This time, the wave functions of the type (iii) and of the type (iv) correspond 
respectively toe = = O, e = = 1/2 and e = = 1/2, e = = 0. After some 
algebra, one can find the corresponding operators H and also the matrices Mi, Mi, M 
are deduced. For these two types (iii) and (iv), we got the three same QES conditions 

CC 

M ■ M I// • 2 , MM - An{n + l) = (30) 

I6p 2 

The condition involving M fixes the difference between the constant A, A, namely 

A — A = M 2 — M 2 . (31) 

It is found that this QES condition is the same for all four types of the wave function. 

A consequence of these results is that the equations (17) admits a double alge- 
braisation. Solutions of the types % and ii exist if the condition (28) is fulfilled and 
solution of the types iii and iv if the condition (30) holds. 
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M 

Figure 1: The critical value of p as a function of the coupling constant M for the type 
i solution and n = 2. The integers label the number of real algebraic eigenvalues 

3.2 Some properties of the spectrum 

We studied the algebraic eigenvalues of the equation (17) for the solution of the type 
i and for n = 1,2. The invariant vector space has dimension 2n keeping into account 
that the condition (5) imposes a constraint on the polynomials. In the case n = 1 
the two algebraic eigenvalues have the form 

E = p 2 + 2 ± ^l-p^l + M) 2 . (32) 

Showing that the algebraic eigenvalues are real only for |p| < 1/|1 + M\. Similar 
features are observed namely in [7], [8]. Namely the eigenvalues come out to be real 
or in complex degenerate pairs. For n = 2 the polynomial giving the four algebraic 
eigenvalues for E is real but rather involved and the solution cannot be expressed in 
a closed form for generic values of M, p. In the case M = 1 and M = 3, however, we 
could find explicit solutions : 

M = 1 , E = 4 + p 2 (2 times) , E = 8 + p 2 ± 8^1 - p 2 

M = 3 , E = 10 + p 2 ±^9- V , E = 2 + p 2 ±2 v /l-4p 2 

The plane M, p is partitionned into regions admitting 4,2 or algebraic eigenvalues. 
The critical values of p c are presented as functions of the parameters M for < M < 3. 
We note in particular that one of the critical line become infinite in the limit M — > 1, 
indicating that in this case there are two real eigenvalues for p > 1. 
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3.3 Hyperbolic case 

The construction of previous section can also be realized for the case where the 
trigonometric functions entering in the potentials are replaced by their elliptic coun- 
terpart. The discussion of the different algebraisations turn out to be the same. Here 
however, we will study in detail the algebraic properties of the operator given by 



i,n-o,h±r -;Uf C(cosh2x L) • C 
\ £>(cosh2:r- 1) + D —£s - (pcosh2:r - iMf j ' 



dx 2 



where p is a free real parameter. One can assume M > M without loosing generality. 
The gauge transformation is performed as follows 



H = exp(— 6cosh2x)H exp(# cosh 2a;), 
H n H12 
if 21 H22 



(34) 



On further substituting z = cosh 2x — 1 and fixing the constant 9 by means of 9 = y , 
the different components of H read 

~\ d 2 , . , d d 9 J , 9 d M — 1 

H n = -4,( 2 + 2)- - 4(* + 1)- - S, PZ - - f - - — Z ) 

+2ip{M- 1) + M 2 , 
H 12 = Cz + C, 
H21 =Dz + D, 

~\ d 2 w \ d „ . d 9 , . , o d M — 1 . 
tf 22 = - 4z ( 2 + 2)^ - 4(* + 1)- - 8 V2 - - <? - 4^- - — ,) 

+2ip(M - 1) + M 2 . (35) 

Decomposing now the operator H according to Eq.(l), we obtain 

-Aip{z 2 f z - Nz) Cz 

Dz -4ip(z 2 £ - Nz) 



^ ^ , (36) 



where we posed N = N = The form of Hq and can be obtained 

easily. Note that (i?i)r2 = and (#1)21 = Dz, so that 5 = 5' — 1 in this case. 
Refering to the above general case and with 

fa z n \ fa lZ n - 



* = ^J + U>' + -- (37) 



we can write the vector Hip according to 



Hip = diag(z n+1 , z n+1 )M 1 



a 



+ \diag(z n ,z n )M 1 (^) +diag(z",^)M ^ 
+ ..., (38) 
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where 



,-Aip(n-N) C 
Ml ~\ D -Mp{n - N) 

f-4ip(n-l-N) C 
1 " V D -4ip(n-l-N) 

Mo =- (to . + «^ + ^). + (^H^ + i)' itp , + % + l)2 ).m 

The three necessary conditions for the operator H to have a finite dimensional 
invariant vector space can then be obtained in a straightforward way, the final results 
read : 

N + N = 2n-1 , 16p 2 (n-N)(n-N) + CD = , ^ = 4 ^( w ~ ff) ( 40 ) 

the equation involving the metric M imposes in turn 

Oft + 4(N - N) (2n + ip)[3 a - Da 2 = (41) 

As a result, assuming a choice of the integer n, we end up with a family of QES 
operators labelled by the parameters N, p, C/D and C. 

Different choices of the non-diagonal interactions H 12 and H 2 \ can be performed 
which lead to similar conditions between the cosmological constants. We will discuss 
these possibilities in the framework of periodic potentials (formulated in terms of 
trigonometric functions), largely discussed in the next section. 



4 PT-symmetric QES equation with polynomial 
potential 

In this section, refering to unidimensional operator studied in [9] we will construct a 
PT-symmetric QES matrix Hamiltonian of the form 

H = --^H 2 + M 4 (x) (42) 

where M^(x) is 2 x 2 PT-symmetric matrix. The above Hamiltonian can be written 
in terms of components and we choose the potentials of the form : 

d 2 

H u = — r^7 - x + iAx3 + Bx 2 + iCx + A 
ax 2 

H X2 = u, 

#21 = U>, 

A 1 .... 

H 22 =-—-x 4 + iAx 3 + Bx 2 + iCx + D (43) 
ax 2 
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In order to reveal the QES property, it is convenient to perform a gauge transformation 
according to 

H = exp(— ax 3 — (3x 2 — r yx)H exp(aa; 3 + (3x 2 + 72;), (44) 
The gauged Hamiltonian then simplifies considerably if 

cl = ~, P = -j, 1 = X -{B~),A = A,B = B. (45) 
leading to the following expression : 

H = — — — 4f3x-^ — 27- 6a[(x 2 - mx) + 9xa 3 ] 

ax 2 ax ax ax 

+ {-2(3 - 7 2 ) + diag(D, D) + uoa + + u)<7_, (46) 

where the constants C, C have been redefined according to C = i(6a(m — 9) + 6a + 
4/3 7 ), C = i(6a(m + 9) + 6a + 4/3 7 ). 

However, in this form, the occurence of an invariant finite dimensional vector space 
of function is not yet manifest in the sense that the operator H doesn't preserve the 
vector space {P m ~e, Pm+oY- 111 order to reveal such a possibility we can apply the 
technique of the first section. Here we will follow [5] and perform a supplementary 

transformation on the operator H with the matrix S = f ® x J . After an algebra, 

we obtain finally the form 

H =S- 1 HS, 

-[-^ + ^-^-t)s +d+ i( b -t) 2 i-<- 

A A d 2 
+2zdiag( J + (n - 2), J + (n)) + diag(-, --) + ua_ - uj\ 2 -^a + (47) 

with J+(n) = x 2 -j- — nx. Here we have set m — n — 1,D — —A + D,9 = 1 and 
fixed the arbitrary parameter A entering in the gauge transformation by means of 
uj = —H\n. 

The Hamiltonian H manifestly preserves the finite dimensional space (P n _2,i 3 n )*. 
Note that A,B,D,Co are free real parameters, n is a non-negative integer and A is a 
free complex parameter. 



5 Recurence relations 

In this section we will express the formulation of the QES solution in terms of re- 
curence relations to the case of PT-symmetic matrix Hamiltonian. We will see that 
the eigenvalue equation Hip = Eip leads to a system of four terms recurence relations. 
The solutions ip are of the form 

if \ i ix3 Ax2 1 »/d A \ \ Pk(E)x k \ ,. Q . 

^ {X) = 6XP( "T - — + 2 {B - ~A )X) ( ZZMEW ) (48) 
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To solve the equation Hip = Eip is equivalent to solve the following equation 

\ E? =0 Qi(E)x l ) fi lE™o QiW ) 1 J 

Then the equation above can be transformed into a fourth-order recurence relation. 
It reads 



where 




(k{k-l) 

-Co (k + 2)(k + l) 



A? 
4 



)+Mv](k-l) 
l-\u + i(B-^-)}(k + l) 



-D-\(B-f) 2 -A(k-2)-j + E Co\ 2 k(k-l) \ 

-D- \(B-f) 2 -Ak + j + Ej 

In the present case, the recurence relations are of fourth order, contrasting with 
other cases studied in the litterature [10, 13] where they are of third order. Setting 
u = uj = the two recurence relations decouple and the corresponding equations (e.g. 
the one for P^) correspond to the scalar PT-invariant and QES quartic oscillator. It 
is also of fourth order; as a consequence, both P and Pi are arbitrary (P fixes the 
normalisation) and the other P^, k > 2 are determined recursively. The construction 
of the QES eigenvalues associated to this system is not as transparant in in the case of 
third-order recurence where a common factor, say P n factorize out of the Pfe's, k > n. 
In the present case, the QES eigenvalues are obtained by solving the system 

P n (P,P) = , P n _ 1 (P,P 1 )=0 (52) 

which is linear in P x . These conditions indeed lead to a truncation of the series for 
ipi(x) defined in (48). Coming back to the full system (i.e. with uj ^ 0, Co ^ 0), it 
is easy to see that Qq 1 Q\ 1 Q2iQz remain arbitrary (Q set the normalisation). The 
QES eigenvalues can be obtained by solving the system 

P n (E,Q 1 ,Q 2 ,Q 3 ) = , P n -i(E,Q u Q 2 ,Q 3 ) = 
Qn+2(E, Qi, Q2, Q3) = , Q n+ ±(E, Q 1 , Q 2 , Q3) = 

(53) 

which turns out to be linear in Qi, Q 2 , Q3. 
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6 Conclusions 



In this paper, we have proposed a set of simple nesessary and sufficient conditions for 
matrix-valued operators of a certain type to preserve a vector space of polynomials 
of fixed degrees. We have seen that the scalar Razhavi potential admits QES matrix 
extensions of several types. We also constructed a QES, matrix- valued PT invariant 
Hamiltonian with polynomial potentials. Finally, by taking this last problem as 
example, we have shown that the coupled differential equations can be transformed 
into a system coupled recurence equations of fourth order. 
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